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ABSTRACT 

TeV photons from blazars at relatively large distances, interacting with the optical-IR cosmic 
background, are efficiently converted into electron-positron pairs. The produced pairs are ex- 
tremely relativistic (Lorentz factors of the order of 10^-10^) and promptly loose their energy 
through inverse Compton scatterings with the photons of the microwave cosmic background, 
producing emission in the GeV band. The spectrum and the flux level of this reprocessed 
emission is critically dependent on the intensity of the intergalactic magnetic field, B, that can 
deflect the pairs diluting the intrinsic emission over a large solid angle. We derive a simple re- 
lation for the reprocessed spectrum expected from a steady source. We apply this treatment to 
the blazar lES 0229H-200, whose intrinsic very hard TeV spectrum is expected to be approxi- 
mately steady. Comparing the predicted reprocessed emission with the upper limits measured 
by the Fermi/Large Area Telescope, we constrain the value of the intergalactic magnetic field 
to be larger than B '^'a x 10^^^ Gauss, depending on the model of extragalactic background 
light. 

Key words: galaxies: active - galaxies: jets - BL Lacertae objects: individual: lES 0229H-200 
- radiation mechanisms: non-thermal - gamma-rays: observations. 



1 INTRODUCTION 

Although cosmic magnetic fields are present at all scales and likely 
play a non negligible role in shaping the Universe, our knowledge 
of them is still rather poor. In particular, magnetic fields at cos- 
mological scales (> 1 Mpc) remain elusive (e.g. Kronberg 2001, 
Neronov & Semikoz 2009). From the theoretical point of view, 
intergalactic space can contain the traces of magnetic fields pro- 
duced in the initial phases of the Universe (e.g., Grasso & Rubin- 
stein 2001) and can also be polluted by magnetic field ejected by 
galaxies and quasars (e.g. Furlanetto & Loeb 2001). Methods based 
on the rotation measure in the radio band constrain the intensity of 
the intergalactic magnetic field (IGMF) below B < 10""-10"'-^ 
G (e.g. Kronberg 2001 and references therein) but these values de- 
pends on the highly uncertain correlation length of the field. 

As early proposed by Plaga (1995), IGMF could be probed 
by using indirect methods based on the effects of the IGMF on the 
electron-positron pairs produced through the interaction of a TeV 
photon from a cosmological source with a low energy photon of the 
optical-IR cosmic background. The produced pairs promptly loose 
their energy through inverse Compton (IC) scattering, in which a 
photon of the cosmic microwave background (CMB) is boosted at 
7-ray (GeV) energies. The existence of a magnetic field, even with 
the tiny intensities suggested by the measures mentioned above. 
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modifies the properties of the 7-ray emission, causing observable 
effects (time delays, e.g. Dai & Lu 2002; formation of extended 7- 
ray halos, Aharonian, Coppi & Voelk 1994, Neronov et al. 2010) 
that can be exploited to infer the value of B (e.g. Neronov & 
Semikoz 2009 and references therein). The reprocessing of multi- 
TeV emission from blazars into GeV-MeV emission has been also 
proposed has a possible contribution to the observed 7-ray back- 
ground (Coppi & Aharonian 1997, Venters 2010). 

In this paper we use a method based on the observed level of 
the GeV emission of an highly absorbed BL Lac, lES 0229+200, 
whose intrinsic TeV spectrum is expected to be almost steady and 
very hard (Aharonian et al. 2007, Tavecchio et al. 2009). Several 
papers discussed a similar method applied to the emission from 
gamma-ray bursts (e.g. Dai & Lu 2002, Razzaque et al. 2004, 
Takahashi et al. 2008). A similar method for blazars have been al- 
ready outlined in Dai et al. (2002) and Murase et al. (2008), but 
these works were focused on the discussion of the variable GeV 
emission, by-product of the reprocessing of the primary TeV emis- 
sion of a rapidly flaring blazar. Here, instead, we present a sim- 
plified treatment of the case of stationary emission of the primary 
TeV source, suitable to describe the case of lES 0229-1-200. Basi- 
cally, our approach is based on the fact that for increasing values of 
the IGMF the level of expected reprocessed emission in the MeV- 
GeV band decreases, since the total amount of the primary flux is 
spread over larger solid angles. The comparison between the ex- 
pected level of the reprocessed flux and the flux at GeV energies 
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Figure 1. Cartoon of the reprocessing of the absorbed TeV radiation. The 
source illuminates the regions inside a cone with semi-aperture dc- TeV 
photons are converted into e+/^ pairs at a typical distance of hundreds of 
Mpc. Pairs cool rapidly through inverse Compton scattering on the CMB. 
For a given intensity of the (perpendicular component of the) extragalac- 
tic magnetic field, pairs with very high energy (upper panel) cool without 
changing their direction (arrows) and thus the reprocessed GeV emission 
from the IC is beamed within the same angle d c- The emission of progres- 
sively lower energy of the pairs is instead spread over larger angles, due 
to the curved trajectoiies of the pairs, in turn due to their longer lifetimes. 
As detailed in the text the resulting overall spectrum (right panels) can be 
approximates by three power laws, each branch due to electrons emitting 
within the oiiginal cone (upper), within a cone larger than the original (mid- 
dle) and almost isotropically (lower). 



observed by the Fermi/Large Area Telescope thus provide a value 
(or a lower limit) on the intensity of the IGMF. 

We use the notation Q — Qx 10^ in cgs units. 

After the completion of this work a paper on same argument 
have been published (Neronov & Vovk 2010). The authors derive a 
value of the magnetic field very close to the value derived here, us- 
ing a full treatment considering also the development of a cascade. 
However, they implicitly assume that the source is emitting TeV 
radiation isotropically. In this case the total flux of the reprocessed 
emission does not depends on the IGMF. The IGMF has the only 
effect to reduce the expected surface brightness of the reprocessed 
emission, lower for larger values of the magnetic field. Instead, in 
our work we assume that the primary high-energy photons from the 
blazar, as commonly assumed, are strongly beamed since they are 
produced in a narrow relativistic jet. In this case the effect of the 
IGMF is to spread the reprocessed emission in solid angles larger 
than the original beaming angle and thus reduce the total observed 
flux, even if the source remains unresolved. 



2 THE REPROCESSED SPECTRUM 

The interaction of a 7-ray with energy E = 1 EtcV TeV with a 
low frequency (optical-IR) photon of the extragalactic background 
light (EBL) produces an electron positron pair with Lorentz factor 
7 ~ E/2me(? — IO^-EtcV (for simplicity and since the redshift 
of lES 0229-1-200 is small, we neglect here and in the following the 
redshift dependence of the photon energies). These pairs, in turn, 
will inverse Compton scatter the photons of the CMB producing 
7-rays of (observed) energy e ~ ^^hvcMB — 2.8 fcTcMB 7^ = 
0.63 ^Tev GeV. In the following we assume that the maximum en- 
ergy of the photons produced through IC scattering, emax, is below 
the threshold of absorption by the EBL (approximately above 400- 
500 GeV for z = 0.14 for most of the EBL models). For higher 
Emax an electromagnetic cascade would develop and a full numeri- 
cal treatment is needed (e.g. Coppi & Aharonian 1997, d' Avezac et 
al. 2007; but see §4 below). 

The typical distance from the blazars at which 7-rays are ab- 
sorbed is of the order of few hundreds Mpc for 1-10 TeV photons 
(e.g. Dermer 2007). On the other hand, once produced, pairs cool 
very rapidly, the cooling time for IC scattering being: 

tcool = . ^, ^"^'^ ^ , ,4 ^ 7 X 10^Se'^(l + ^O'^S (I) 

47(JTtVcMB,o(l + Zr)* 

translating into a distance rfcooi = 2 x 10^*7^^(1 + Zr)~* cm. 
Here (7cmb,o is the energy density of the cosmic background radi- 
ation evaluated at 2; = 0, and Zi is the redshift of the reprocessing 
region (not necessarily equal to the redshift of the primary source). 

Therefore, at least for the pairs with the largest energies, we 
can approximate the conversion and emission zone as a thin slice 
(see also Ichiki et al. 2008). The electrons will have a completely 
cooled distribution. A*' (7) oc 7"^, and the corresponding repro- 
cessed emission will be characterised by a spectrum f (e) oc e^"'^. 
This power law extends up to an energy emax determined by the 
maximum Lorentz factor of the produced pairs, in turn depending 
on -Bmax, the maximum energy of the 7-rays emitted by the blazar, 
==;63(£;max/10TeV)2GeV 

The geometry we adopt is sketched in Fig. [T] Since the emis- 
sion from the blazar is strongly collimated into a cone with semi- 
aperture ~ 1 /r, where P is the bulk Lorentz factor of the jet, 
we can approximately assume that the source (located at the apex 
of the cone) uniformly illuminates the gray conical region in Fig.[T] 
Photons are converted into pairs that produce the GeV emission at 
the base of the cone. This geometry is different than that assumed in 
other works (e.g. Neronov & Vovk 2010) that assume that the pri- 
mary TeV emission is isotropic. In this case the total flux recorded 
from the source does not depend on the intensity of the IGMF: only 
the surface brightness depends on it, since the effect of a larger 
IGMF is to broaden the (energy dependent) size of the reprocessed 
"halo". Instead, in our scenario, the extension on the sky of the re- 
processed emission cannot be larger than 0^ « 9cd^/{D — d^) 
where d-, is the distance (from the source) where the 7-7 optical 
depth is unity, and D is the Earth-source distance. 

The reprocessed emission contains all the absorbed flux. In 
absence of any magnetic field, B = 0, the pairs are produced and 
emit along the direction of the primary 7-ray (Fig.[T] upper panel) 
and the observer would measure a total reprocessed GeV flux equal 
to the absorbed TeV one (although in a different energy range). The 
condition that the absorbed flux is completely reprocessed as IC 
radiation is then: 



Fint{E)[l- 
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where Fint(-E) is the source intrinsic spectrum (extending from 
-Emin to -Bmax) and T (E) IS the energy-dependent optical depth 
of 7~rays. Further assuming: 

-0.5 



F{e) = k 



e < En 



(3) 



Eq. |2]can be used to derive the normalisation k = Fabs/(2emax) 
(for Eynin < -Bmax, implying also emin < tmax)- We remark that, 
if the intrinsic blazar spectrum is harder than E^^, the value of 
the absorbed flux Fabs (and, in turn, the level of the reprocessed 
radiation) depends solely on the maximum energy of the intrinsic 
emitted spectrum, -Emax- 

The opposite case is that of a magnetic field so large that all the 
produced pairs are promptly isotropized in a time shorter than tcooi 
(bottom panel in Fig. [T](. The total flux recorded by the observer 
will be therefore strongly diluted. Note that in the isotropic case 
the observed reprocessed emission contains the contribution of the 
two jets, that pointing toward us and that pointing far from our line 
of sight. As above we can write: 



f(e)de = 2Fab, 



4-K 



(4) 



where Qc — t^d't is the solid angle in which the primary radiation is 
coUimated. If the emitting source in moving with V in one direction, 
then Q,c ~ tt/F^. The factor 2 takes into account the contribution 
of both jets. The reprocessed spectrum is still given by Eq.[3]but 
with a normalisation k — -Fabs^^c/(47remax). 

More complex is the intermediate case, for magnetic field such 
that electrons with different Lorentz factors are deviated by dif- 
ferent angles. In this case, the reprocessed emission is beamed 
into different solid angles Q~f — 27r(l — cos 9-,) at different 
energies. The angle 9^ can be estimated assuming it is the an- 
gle by which the electron velocity vector changes in the cooling 
length ctcooi: 9j ~ ctcooi/rL = 1.17B-i5(l + Zr)~*7ff^ where 
tl = jmeC^ / (eB) ~ 2 x 10'^*76-BZi5 cm is the Larmor radius 
of the electron. 

Since 9j oc 7~^, we can derive the slope of the emitted spec- 
trum in this regime: 



F(e) oc 7 oc e ' 



(5) 



where we used A''(7) oc 7~^, 7 oc 7^ (IC scattering), and e oc 7^. 
Therefore the reprocessed spectrum in this region is described by a 
very hard power law (middle panel in Fig.[T}. 

In real applications, there will be a combination of all the three 
cases described above: for a sufficiently low energy, e < eiso the 
Lorentz factor of the corresponding emitting electron is sufficiently 
small to be isotropized by the magnetic field B in its cooling time 
and the spectrum around e will be thus given by Eq.|4l For energies 
above eiso but lower than an energy to electrons curves substan- 
tially under the action of the magnetic field, but the emission is not 
isotropic and the spectrum is described by the hard power law e^^^. 
At energies above (the energy at which approximately 9^ ~ 9c), 
the reprocessed radiation is insensitive to the magnetic field, since 
the curvature of the electrons during their emission is (much) less 
than the beaming cone of the intrinsic blazar emission. The total 
spectrum will be described by three power laws and it can be well 
approximated by the following relation: 



F(e) = k(^\ 

\ ^max / 



e < En 



(6) 



The normalization k is found in this case noting that Eq. |6]has to 



reduce to Eq.[3]for the case B — 0. Inserting then Q,-, — above 
and comparing the resulting expression with Eq. [3] we derive the 
value of the normalization: k — r2c-F'abs/2emax. 

For simplicity, in Eq. |6] we do not include the contribution 
from the receding jet, visible for energies at which the pairs are 
almost completely isotropised. This contribution would produce a 
doubling of the flux in the isotropic case, in a region of the spectrum 
not interesting for the case we are discussing here. 



3 APPLICATION TO lES 0229+200 
3.1 LAT data 

lES 0229-1-200 is not present in the list of bright AGN from the first 
three months of Fermi Large Area Telescope (LAT) observations 
(Abdo et al. 2009). Therefore, we searched for detections in the 
publicly available datfl We selected the photons of class 3 (dif- 
fuse) with energy in the range 0.1-100 GeV collected from 2008 
August 4 (MID 54682) to 2010 March 31 (MJD 55286), for a total 
of 19 months of elapsed time. These data were analyzed by us- 
ing Science Tools 9.15.2, which includes the Galactic diffuse and 
isotropic background and the Instrument Response Function IRE 
P6 V3 DIFFUSE. After a selection of the useful events and good- 
time intervals, specifically taking into account a zenith angle < 105 
deg to avoid the Earth albedo and the photons within the region of 
interest (ROI) with radius of 10 deg from the source radio position, 
we prepared a count map of the ROI. Then, we searched on the 
map for other possible contaminating sources inside the ROI, be- 
cause we have to take into account them in the model. We model all 
the sources as single power laws, with flux and photon index free 
to vary. The following steps are to calculate the live-time, the ex- 
posure map and the diffuse response. With all these information at 
hands, we performed an analysis by using an unbinned likelihood 
algorithm and calculated the corresponding test statistic (TS, see 
Mattox et al. 1996 for a definition of TS; here it is worth noting 
that the significance sigma roughly equal to VtS ). None of the 
energy bins shows a significance larger than TS — 25, therefore 
only upper limits can be derived (although for the 10-30 GeV bin 
there is a hint of detection at the VTS ~ 4(t level). Following 
Abdo et al. (2010a), upper limits have been calculated looking at 
the flux satisfying 2Alog(I/) = 4 when increasing the flux from 
the value minimizing the likelihood L. 

The measured flux and photon index are the average over 
the 19 months of analysed data. Since the SED studied here are 
built with non simultaneous data and the lES 0229-1-200 is a weak 
gamma-ray source, we have not studied the variability of the de- 
tected sources and used only the averages. The results are sum- 
marised in Table I. The latest available calibration of LAT indi- 
cates the presence of systematic errors, which have to be added to 
the quoted statistical errors in Table 1, and have value of 10% at 
100 MeV, 5% at 500 MeV, and 20% at 10 GeV (Rando et al. 2009). 



3.2 Results 

Fig. 2 shows the high-energy spectral energy distribution of lES 
0229-1-200 including the TeV data from H.E.S.S. (red) and the 
same points corrected for the absorption (green) using the Low 
SFR model of Kneiske et al. (2004). We assume that the intrin- 
sic spectrum is well represented (dashed line) by a hard power law. 
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Figure 2. Spectral Energy Distribution of lES 0229+200 in tlie liigli energy band and tlie expected reprocessed GeV emission. Red points sliow the observed 
H.E.S.S. spectrum (Aliaronian et al. 2007) and the green points the points after the connection for the absorption with the LowSFR EBL model of Kneiske et al. 
(2004). See Tavecchio et al. (2009) for details. The short-dashed black line is an approximation of the intrinsic spectrum, modelled as a hai'd power law. The 
long-dashed line is the con'esponding absorbed spectrum following Kneiske et al. (2004). The shaded gray region between these two lines shows the absorbed 
flux. Note that the total amount of absorbed power, that is the relevant quantity for the estimate of the level of the reprocessed GeV emission, is only slightly 
sensitive on the assumed spectral shape, since it is dominated by the flux at the largest energies. The coloured lines report the approximation of the expected 
reprocessed spectrum for different values of the IGMF, B = 10^^'', 10~^*, 10"^'' G and two different values of the initial collimation angle (9c = 0.1 
rad and 0.05 rad, solid and dashed-dotted lines, respectively), determining the intrinsic beaming of the primary radiation. For comparison, we also report the 
curves (black dotted lines) corresponding to the two extreme cases of (upper) B = and (lower) completely isotropy of the reprocessed emission (extremely 
large B). Black points show the Fermi/LPiT upper hmit to the flux in the 100 MeV-1 GeV, 1-10 GeV and 10-30 GeV bands obtained through the spectral 
parameters derived with the standard analysis (see text for more details). The solid black line shows the lowest possible reprocessed spectrum consistent with 
the upper hmits, coiTesponding to a magnetic field of B = 5 X 10^^^ G (for 8c = 0.1). 
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3 


1.55 
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17.65 



Table 1. Results of the analysis of the LAT data. (1) Energy band, in GeV. 
(2) Estimated upper limit (ph cm"^ s" i). (3): number of photons in the 
band. (4): value of TS. See text for details. 

Fe oc E^^'^ (see Tavecchio et al. 2009 for the justification of this 
choice). The black dotted line is the corresponding absorbed spec- 
trum using the Kneiske et al. (2004) model. The area in gray show 
the flux absorbed and available for reprocessing. As long as the in- 
trinsic spectrum is hard, the amount of absorbed energy depends 
only on the intrinsic luminosity of the highest energy bin iSmax- 
The most conservative limit on the IGMF corresponds to the low- 
est amount of reprocessed radiation that in turn corresponds to the 
smallest intrinsic luminosity. To this aim we use the EBL model 
providing the lowest opacity around 10 TeV. 

We report the LAT upper limit in the 0.1-1 GeV, 1-10 GeV 



and 10-30 GeV band (Table 1). The solid and dot-dashed lines 
report the expected reprocessed emission assuming three different 
values of the IGMF and two different beaming angles (0.05 and 0. 1 
rad, corresponding to bulk Lorentz factor F = 20 and 10, respec- 
tively) for the intrinsic blazar emission. The black line is calculated 
for the minimum value of the magnetic field consistent with the 
upper limit, B = 5 x 10-^^ G. Note that, due to the very hard 
reprocessed spectrum, the most stringent upper limit is that at the 
highest energies, 10-30 GeV. Beaming angles smaller than those 
assume here (corresponding to larger bulk Lorentz factors of the 
jet) would result in lower values for the upper limit on B, see Eq. 

m 

We remark that, unlike the methods based on the estimate of 
the rotation measure in the radio band (e.g. Kronberg 2001), with 
which it is possible to derive upper limits to the IGMF, this methods 
allows us to put a lower limit on B. If the hint of detection in the 
highest energy bin is real, we have two possibilities: either it is 
the reprocessed radiation, and in this case it gives a measure of B, 
or it is still primary emission from the blazar (even if belonging 
to a component different than that observed at TeV energies, e.g. 
Tavecchio et al. 2009), and in this case our limit would still hold. 
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4 DISCUSSION 

The lower limit on the value of the magnetic field derived here, 
can be considered one of the most stringent value 
ever derived for the IGMF. The value is mainly constrained by the 
LAT upper limit above 10 GeV, in which the source is tentatively 
detected at the 4a level. We remark that, since we derive a lower 
limit on B, even if this FermifLKT measure is considered as an 
upper limit the conclusion does not change. 

The main assumptions adopted in this paper are: (i): the 
amount of reprocessed energy is derived from the level of the ob- 
served spectrum measured by H.E.S.S. (Aharonian et al. 2007); (ii) 
we use the Low SFR Rneiske et al. (2004) model to calculate the 
optical depth; (iii) we assume that the maximum energy of the in- 
trinsic absorption is 11 TeV. 

Assumption (i) is supported by the extremely small variability 
observed by H.E.S.S. during the observations in 2005-2006. Of 
course, any variations of the intrinsic emission of the blazars, will 
be reflected into the reprocessed emission. Ideally, simultaneous 
TeV and GeV observations would be required to take into account 
variability and exploit the information carried by the delay between 
the intrinsic TeV and the reprocessed GeV emission (see Dai et al. 
2002, Murase et al. 2008 for discussions). However, this approach 
is prevented in the specific case of lES 0229-1-200 by the very small 
GeV flux, still at the detection limit of LAT after 18 months of 
observations. Also in other TeV BL Lacs more bright in the TeV 
band, integration over few days are necessary to detect the GeV 
emission (e.g. Abdo et al. 2010b). Moreover, in these objects, the 
intrinsic GeV emission is much more luminous than the expected 
reprocessed radiation. These features make difficult to reveal the 
reprocessed radiation even during bright flares. 

The use of a specific model for the absorption is unavoid- 
able. However, several of recent EBL models converge at the wave- 
lengths above 10 fim, the ones determining the opacity for 7-rays 
of energy below 10 TeV. In particular, the recent models of Frances- 
chini et al. (2008), Gilmore et al. (2009), Finke et al. (2010) agree 
well (see e.g. the discussion of Finke et al. 2010) and are consistent 
with the low level of the EBL suggested by recent observations of 
Cherenkov telescopes (e.g. Aharonian et al. 2006, Mazin & Raue 
2007; see also Kneiske & Dole 2010). The model adopted here 
(LowSFR of Kneiske et al. 2004) provides an optical depth similar 
to that of all the other updated models up to energies of 4—5 TeV. 
For larger energies the predicted optical depth is lower than that of 
the other models (see also Tavecchio et al. 2009). We again stress 
that this implies to minimize the luminosity of the reprocessed ra- 
diation, and hence the derived lower limit of iJ is a conservative 
value. 

The assumed maximum energy of the intrinsic emission, point 
(iii), is somewhat critical, since it determines the total amount of 
energy reprocessed and then re-emitted in the GeV band, and the 
possible development of cascades. The first point is clear: since 
the intrinsic spectrum is rather hard, the total luminosity of re- 
processed radiation (and its maximum energy) depends on -Bmax. 
Under the assumption that the spectrum is proportional to E^^^, 
the normalization of the reprocessed emission is (see Eq. |6ll fc oc 
Smax/emax oc -Emax. Therefore all the curves reported in Fig[2] 
will shift upward by this factor. We consider unlikely that the spec- 
trum extends unbroken at energies well above few tens of TeV and 
therefore the correction factor cannot be much larger than a factor 
of a few. In any case we remark again that, since a larger -Emax 
would have the effect to increase the value of the derived lower 
limit (a larger B would be required to downshift the curves), our 



limit can be considered also in this respect, again, as a conservative 
value. The same argument applies when considering the possible 
role of electromagnetic cascades. For values of -Bmax larger than 
~ 30 TeV a sizeable amount of the reprocessed flux would be emit- 
ted above ~ 500 GeV (the energy above which the optical depth for 
pair conversion exceeds 1) and then reabsorbed, in turn producing 
new pairs and thus initiating a cascade. In this case, apart the differ- 
ent resulting spectrum, the total flux of the reprocessed component 
would be larger than what simply estimated here: therefore, also in 
this case, our value of B can considered a conservative value. 

Finally, we note that in our derivation we are implicitly as- 
suming that the magnetic field is oriented perpendicularly to the 
direction of the relativistic pairs. In reality the IGMF is probably 
randomly oriented, maintaining its coherence in domains with a 
size ~1 Mpc (e.g. Neronov & Semikoz 2009). Including the ge- 
ometry of the field in our calculations would result in a limit to 
the IGMF slightly larger than that derived here. Again, our inferred 
value of B is then a conservative one. 
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